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Abstract The conventionally separated treatments for 
strangelets and strange stars are now unified with 
a more comprehensive theoretical description for ob- 
jects ranging from strangelets to strange stars. Af- 
ter constraining the model parameter according to 
the Witten-Bodmer hypothesis and observational mass- 
radius probability distribution of pulsars, we investigate 
the properties of this kind of objects. It is found that 
the energy per baryon decreases monotonically with in- 
creasing baryon number and reaches its minimum at the 
maximum baryon number, corresponding to the most 
massive strange star. Due to the quark depletion, an 
electric potential well is formed on the surface of the 
quark part. For a rotational bare strange star, a mag- 
netic field with the typical strength in pulsars is gener- 
ated. 
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1 Introduction 


It was pointed out long ago that strange quark mat- 
ter (SQM) might be the ground state of strongly in- 
teracting matter, which is nowadays called the Witten- 
Bodmer hypothesis {1|/2|. If true, there should exist sta- 
ble objects of SQM with the baryon number A ranging 
from a few to ~10°". Customarily, small SQM nuggets 
with A < 10’ are often referred to as strangelets [Ho], 
or slets mp while stars consisting of SQM are called 
strange (quark) stars uže], being possible candidates 
for pulsars. 

Lumps of SQM are expected to be produced in the 
collision of binary compact stars containing SQM 


[23]. Further collisions among those lumps may create 


slets, nuclearites p423], meteorlike compact ultradense 
objects etc., and some of them may eventually make 
their way to our Earth 27]. Due to the special char- 
acteristics of these objects such as the lower charge- 
to-mass ratio (28][29], the larger mass [B0], the highly 
ionizing tracks in the interstellar hydrogen cloud (e.g., 
pulsar scintillations) , and the characteristic gamma 
rays through heavy ion activation [82], there are pos- 
sibilities to observe them. However, despite decades of 
efforts, no compelling evidence for the existence of sta- 
ble SQM is found (for reviews, see, e.g., Refs. [83][84]). 

This case is due to the extreme complexity of 
an SQM system which involves all the fundamental 
interactions, i.e., the strong, weak, electromagnetic, and 
gravitational interactions. In the conventional theoret- 
ical treatments, significantly different simplifications 
were adopted for slets and strange stars. For a slet, 
electrons were ignored since the Compton wavelength 
is much larger than the size of the quark part [85], and 
quarks were assumed to be uniformly distributed. For 
strange stars, gravity has to be considered. The nor- 
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mal way is to first get an equation of state of SQM by 
assuming the local charge neutrality, and then obtain 
the mass-radius (M-R) relation by solving the Tolman- 
Oppenheimer- Volkov equations. 

However, according to recent studies, effects such 
as the charge screening, electron-positron pair creation, 
and nonzero charge densities in strange stars have im- 
portant implications on the properties of SQM. For ex- 
ample, taking into account the electrostatic effects, Al- 
ford et al. found that, for a small enough surface 
tension, large slets are unstable to fragmentation and 
strange star surfaces fragment into a crystalline crust 
made of slets and electrons. For quark-hadron phase 
transition, the finite-size effect turns out to be very 
important [87H40]. It was shown that the geometrical 
structures may be destabilized by the charge screen- 
ing effect [41]. Due to the electron-positron pair cre- 
ation, an upper bound on the net charge of slets or 
strange stars was found [42]. The local charge neutral- 
ity in compact stars is also in question [a3]. In the case 
of a neutron star, an overcritical electric field was found 
in the transitional region from the core to the crust (44). 
For a bare strange star, an electric dipole layer may be 
formed on the surface and results in an electric field 
of ~10!17719 V/cm {12]. Due to the presence of a criti- 
cal electric field, the electron-positron production may 
be induced and results in some astrophysical observ- 
ables [45]. The mass and radius of a strange star are 
increased by ~15% and ~5%, respectively, if the star 
possesses a net charge on the surface [46]. 

Meanwhile, the possibility of pulsars being strange 
stars may give us an insight into the properties of SQM. 
Up till now, around 2,500 pulsars have been observed 
and among them about 70 pulsars’ masses were mea- 
sured (47)[48}. At the same time, more than 10 pulsars 
provide us the M-R probability distributions with pho- 
tospheric radius expansion bursts as well as quiescent 
low-mass X-ray binaries (49}{53}. If SQM is absolutely 
stable, those pulsars may be strange stars 54], then the 
properties of SQM can be constrained with the M-R re- 
lations. 

In the present paper, we study the SQM system 
ranging from slets to strange stars in a unified descrip- 
tion. After constraining the only model parameter, the 
bag constant B, according to the Witten-Bodmer hy- 
pothesis and the observational M-R probability distri- 
bution of pulsars, it is found that the ratio of charge to 
baryon number of a slet is different from previous find- 
ings, while the size is significantly smaller than that 
of a nucleus with the same mass number. In addition, 
rotation of a bare strange star generates a strong mag- 
netic field with the typical strength in pulsars. This 
paper is organized as follows. In Section the uni- 


fied description for SQM objects is presented, where 
the effects of gravity and electrostatic interactions are 
treated on the macroscopic scale while the strong and 
weak interactions are considered locally. Based on this 
description, the properties of SQM objects ranging from 
slets to strange stars are investigated in Section |3| A 
summary is presented in Section [4] 


2 Theoretical framework 


The internal structure of a spherically symmetric, 
charged, and static object should fulfill the thermody- 
namic equilibrium condition, which can be obtained by 
minimizing the energy of the system for given total par- 
ticle number and entropy. We consider the gravity and 
electrostatic interactions on the macroscopic scale. The 
metric for the SQM sphere reads 


ds? = e’dt? — edr? — r?(d6? + sin? 6d¢7), (1) 


where r, 0, and ¢ are the standard spherical coordinates 
with the metric elements satisfying 


2 
e*=1- 2E Mh (2) 
dv 2G 3 a)? 
r7 arr (r- S) + m : (3) 


Here we use the natural system of units, with G and 
a being the gravitational and fine-structure constants. 
The total mass, particle number, and entropy are ob- 
tained with 


Mr) = | arr” (B+ 0@? mr) ar, (4 
0 

Ni(r) = A 4rni(r eò ?r?dr', (5) 
0 

S(r) = i An S(r')e/? 7dr’. È 
0 


Then the total charge is given by Q(r) = 30; aNi(r) 
with qu = 2/3, da = ds = —1/3, and qe = —1. Based on 
the Thomas-Fermi approximation, the pressure P(r), 
energy density E(r), particle number density n; (1), and 
entropy density S(r) are given locally by incorporating 
both the strong and weak interactions. 

By minimizing the mass M = M,(oo) with respect 
to the particle distribution N;(r) and entropy distribu- 
tion S(r) at the fixed total particle number N;(oo) and 
entropy S(oo), we immediately have 


dui  Q nya Mai dv 

dr 2 fio” / 2 dr’ g 
dT T dv 

dr 2dr’ (8) 


with u(r) and T(r) being the chemical potential and 
temperature. 
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For the local properties of SQM, we adopt the bag 
model and consider only zero temperature, where the 
thermodynamic potential density is given by 


Ar) = Qo({ualr)}, tma}) + B (9) 


in the ideal Fermi-gas approximation. To reach the low- 
est energy, SQM undergoes weak reactions and reach 
the chemical equilibrium 


Hu + He = Hd = Hs. (10) 


Then the internal structure of an SQM sphere can be 
determined by solving Eq. (7). Since electrons are not 
confined by the vacuum pressure, an atom-like structure 
of the SQM system is formed, i.e., a positively charged 
SQM core with a cloud of electrons surrounding it. 
The quark-vacuum interface on the surface of the 
SQM core needs to be treated with special care. We 
consider the number of depleted quarks on the interface 
by adopting the multiple reflection expansion (MRE) 


method 


dNs™f  2g;R 


dp; 37 


mR + 1)arctan (=) (11) 


Pi 


givi R 
MinT ( 
Here p; is the momentum of quark flavor i (i = u, d,s). 
Its upper bound corresponds to the Fermi momentum 
on the surface n;(R) with R being the radius of the 
SQM core. Note that Eq. (11) only gives the average 
number of depleted quarks, while for smaller systems 
shell corrections may be important 55]. Then the en- 
ergy contribution and pressure are given by 


vi(R) danse 
=f ppm ap, (12) 
0 dp; 
and 


psut — 5 dE 


Fasurf 
E; 


(13) 


dV | gust 


Under the influence of gravity, the energy contribution 
to the mass is Mf = yo, Beufen(2)/2, 

The quark-vacuum interface is obtained when the 
pressure of quarks is in balance with the vacuum pres- 
sure, i.e., 


P(R) — P.(R) + P™f =o. (14) 


Then we have the total quark number N, = N(R) + 
ee mass M = M,(oo) + Mf, and charge of the 
core Q(R) = X; qi [Ni(R) + NE"). 


3 Results and discussions 
For a given core radius R, the structure of an SQM 


sphere is determined by solving the differential equa- 
tion (7) under the boundary conditions, i.e., M,(0) = 0, 


5.0E-06 
2.5E-06 
0.0E+00 


Fig. 1 (Color online) The density profiles (in fm~?) of an 
SQM sphere with the core radius R = 1000 fm, where the 
bag constant is taken as B1/4 = 152 MeV. Here no, = 
XL; qini corresponds to the local charge density, while np = 
Deas nq/3 is the local baryon number density. 


Q(0) = Q(œ) = 0, and Eq. (14). To illustrate our 
results, we present a colored contour plot in Fig. [1] for 
an SQM sphere with the core radius R = 1000 fm. 
A rich charge profile on the surface is found. At the 
r > R region, there exists an electron cloud which 
neutralizes the positively charged core and expands by 
~1 A. Respectively, the SQM sphere contains 1.37 x 10° 
u-quarks, 1.48 x 10° d-quarks, 1.26 x 10° s-quarks, 
and 2.74 x 10° electrons, which gives the total core 
charge Q(R) = 2.66 x 10° and mass M = 1.22 x 10! 
MeV. Due to the quark depletion on the core sur- 
face (Nswf = —3.34 x 104, Ns“? = —7.06 x 104, and 
Nst — —1.03 x 10°), there exists a surface charge 
Qf =, qi Ns"! = 3.43 x 10°, while the correspond- 
ing mass modification is MS"! = —2.15 x 108 MeV. 


Note that the current masses of quarks and the elec- 
tron mass are taken as m,, = 2.3 MeV, mg = 4.8 MeV, 
ms = 95 MeV, and me = 0.511 MeV [56], leaving only 
the bag constant B undetermined. For SQM to stably 
exist at zero external pressure, the bag constant should 
meet the requirement of the Witten-Bodmer hypoth- 
esis, which gives 144.37 < B'/4 < 159.26 MeV. We 
take the three typical values B!/4 = 145, 152, and 159 
MeV. It is worth mentioning that if B exceeds the up- 
per bound, SQM is unstable and may only exist at the 
core of a compact star, i.e., hybrid star (13157. 


A full calculation from slets to strange stars is car- 
ried out. In Fig. the energy per baryon is given, 
which decreases monotonously with increasing radius 
or baryon number. When the energy per baryon reaches 
930 MeV, the minimum baryon numbers for absolutely 
stable slets are determined, i.e., Amin = 24, 80, and 
394849, which increases dramatically as B approaches 
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Fig. 2 (Color online) The energy per baryon from slets to 
strange stars. In the inset, the obtained M-R relations are 
compared with the normalized M-R probability distribu- 
tion compiled from the observed values of eight pulsars 
Fig. 10]. 


to its upper limit. In the region with 200 fm < RS 
1 km, the variation of the energy per baryon is infinites- 
imal, and the radius is related to the baryon number 
by R = ro A13 with ro = 0.944, 0.901, and 0.862 fm. 
When R = 1 km, gravity starts to reduce the energy 
per baryon and a minimum value is obtained, corre- 
sponding to the maximum mass and baryon number of 
strange stars. The obtained M-R relations, as shown 
in the inset of Fig. |2| are in good agreement with the 
M-R probability distribution obtained by averaging the 
observed values of eight pulsars Fig. 10]. 

The charge-to-mass ratio in Fig. Bla) is defined as 
fz = Q(R)/A ignoring the electrons surrounding the 
SQM core. The obtained values are larger than previ- 
ous results [3}[58}. For smaller B, the SQM core car- 
ries slightly more charge. Note that when R < 13 fm, 
SQM within slets is positively charged. However, for 
R= 17 fm, as shown in Fig. [1] the SQM carries negative 
charge to compensate the positive surface charge. Then 
an electric potential well for negatively charged parti- 
cles is formed due to the quark depletion on the quark- 
vacuum interface. These particles may be trapped in the 
potential well and give a distinct photon spectrum when 
excited, which have significant implications for the ex- 
perimental searches of SQM. As indicated in Fig. [3{b), 
when R = 10° fm, charges are mostly located on the 
core surface and can be described by a constant surface 
charge density o with Q(R) = oR?. It is found that 
o = 0.0144, 0.0135, and 0.0128 fm~? for B1/4 = 145, 
152, and 159 MeV, which are much larger than the 
upper bound 7 x 107° fm? considering the electron- 
positron pair creation [42]. 

When the SQM core rotates, a magnetic field may 
be generated. For strong enough field strength, the elec- 


10° 10' 10° 10° 10° 10° 
Radius (fm) 


Fig. 3 (Color online) (a): The charge-to-mass ratio of the 
SQM core of slets, which is compared with previous findings 
indicated by the dotted line Fig. 4] and dash-dot-dotted 
line Fig. 2]; (b): The surface charge density of the SQM 
core as functions of radius. 
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Fig. 4 (Color online) The magnetic field of rotating strange 
stars estimated with Eq. (15. The dots depict the inferred 
magnetic fields of pulsars obtained from ATNF Pulsar Cata- 
logue [47|[48}, where the binary pulsars are excluded since 
their magnetic fields may be dampened by mass _ accre- 


tion [59]. 
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